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SECOND SUPPLEMENT TO AN ESSAY ON THE 
THEORY OF SYSTEMS OF RAYS 


Read October 25, 1830. 
[Transactions of the Royal Irish Academy, vol. 16, part 2 (1831), pp. 93-125.] 


INTRODUCTION 


The present Supplement contains the integration of some partial differential equations, to which 
I have been conducted by the view of mathematical optics, proposed in my former memoirs. 
According to that view, the geometrical properties of an optical system of rays may be deduced 
by analytic methods, from the form of one characteristic function; of which the partial differen- 
tial coefficients of the first order, taken with respect to the three rectangular coordinates of any 
proposed point of the system, are, in the case of ordinary light, equal to the index of refraction 
of the medium, multiplied by the cosines of the angles which the ray passing through the point 
makes with the axes of coordinates: and as these cosines are connected by the known relation 
that the sum of their squares is unity, there results a corresponding connexion between the 
partial differential coefficients to which they are proportional. This connexion is expressed by an 
equation which it is interesting to study and to integrate, because it contains a general property 
of ordinary systems of rays, and because its integral is a general form for the characteristic func- 
tion of such a system. The integral which I have given in the present memoir, is deduced from 
equations assigned in my former Supplement; an elimination which had been before supposed, 
being now effected, by the theorems which Laplace has established in the second Book of the 
Mécanique Céleste, for the development of functions into series. The development thus obtained, 
proceeds according to the ascending powers of the perpendicular distances of a variable point from 
the tangent planes of the two rectangular developable pencils which pass through an assumed 
ray of the system, and according to the descending powers of the distances of the projection of the 
variable point upon the assumed ray, from the points in which that ray touches the two caustic 
surfaces. In the case of rays contained in one plane, or symmetric about one axis, the partial differ- 
ential equation takes simpler forms, of which I have assigned the integrals, and have given an 
example of their optical use, by briefly shewing their connexion with the longitudinal aberrations 
of curvature. I hope, in a future memoir, to point out other methods of integrating the general 
equation for the characteristic function of ordinary systems of rays, and other applications of the 


resulting expressions, to the solution of optical problems.* 
WILLIAM R. HAMILTON. 


Observatory, 
October, 1830. 


* [See footnote to p. 165.] 
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SECOND SUPPLEMENT 


Statement and Integration of the Partial Differential Equation, which deter- 
mines the Characteristic Function of Ordinary Systems of Rays, produced by 
any Number of successive Reflexions or Refractions. 


1. Suppose that rays of a given colour diverge from a given luminous origin, and undergo 
any number of successive changes of direction, according to the known laws of ordinary reflexion 
and refraction, at surfaces having any given shapes and positions, and enclosing media of any 
given refractive indices. Let a, 8, y, be the cosines of the angles which the direction of a final 
ray makes with three rectangular axes, and let æ, y, z, be the three rectangular coordinates, 
referred to the same axes, of a point upon this final ray; then a, , y, will in general be functions 
of æ, y, z, such that if u denote the refractive index of the final medium, for rays of the given 
colour, the expression 

p (ada + Bdy + ydz) 
is equal to the differential of a certain function V, of which I have shewn the existence and the 
meaning in former memoirs, and which I have called the characteristic function of the final 
system. The design of the present Supplement, is to point out some new properties and uses of 
this function, resulting from the partial differential equation 


AVAE (dV? (dV, 
(ae) + (aq) +a) = Sa 
which we obtain by eliminating the three cosines a, 8, y, between the three equations 
dV dV dV 
Ge ag =p, a= PY (B) 
by the help of the known relation 
e+ e+=. 
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2. The equation (A) is a particular case of a more general differential equation, for all optical 
systems of rays, ordinary or extraordinary, obtained by eliminating the same three cosines, 
a, B, y, by the same known relation between the three following equations, assigned in my former 


memoirs, * 


avo aV w Vd, 
dx 6a’ dy 458’ dz dy’ 


in which V is the characteristic function of the system, and v is a homogeneous function of 
a, B, y, of the first dimension, representing the velocity of the light, estimated on the hypothesis 
of emission, and differentiated as if a, 8, y, were three independent variables. And the integral 
of (A) is a particular case of a more general integral, extending to all optical systems of straight 
rays, and consisting of the following combination of equations, assigned in my former Supplement :* 
bv du bv 
W+ Veart+y 58+? 35° 
sW_, oy oe 
Ba” Sa TY Sasb * Sasy’ 
SW, oy, oy oy 
Se "Sash 4 Spt * * Biy’ 
SW _ æ Sw + Sy +- Z Sv ® 
by" baby V BBB" * By?’ 
between which the three quantities æ, 2, y, are to be eliminated; W being an arbitrary but 
homogeneous function of these three quantities, of the dimension zero; and the partial differential 
coefficients in which the sign ô occurs, being formed by differentiating the homogeneous functions 
W, v, as if a, 8, y, were three independent variables, In applying these general results to ordinary 


systems of rays, we are to put 
| v= p (a+ B+); 


Sa ov’ 8B ov’ by Dv 
2 2 2 
Be Botte, anG- paha- 


or, (making after the differentiations a? + 8? +% = 1,) 


bv ov dv 
e gaas gaa 5 = oy 


oy = p(l at), Samale) pane (l— 9) 
2y 


* [See pp. 111, 112.] 
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and therefore, 
W + V = p (ax + By + y2), 
ô wW 
ga Me PA (aw + By + y2), 


C 
S = uy — pẹ (ao + By +72), oa 


ò wW 
Sy aa PAET 


This system of equations (C) is one form for the integral * of the partial differential equation 
(A); the quantities a, 8, y, being supposed to be eliminated, and W being an arbitrary function 
of these quantities, of the kind already mentioned. 


Transformation and Development of the Integral. 


3. The system of equations (C) may be transformed into the following : + 


dU _ aU Sine a) A 

pe So HY = JB’ Fai See? aes (D) 
in which JU is a function of the three independent variables a, 8, z, obtained from the function W 
by putting 

U= W— pyz, (E) 

and by considering y as a function of a, 8. Let us now proceed to eliminate œ, 8, between the 
three equations (D), by the theorems which Laplace has given in the second Book of the 
Mécanique Céleste, for the development of functions into series. 


This elimination may be simplified by a proper choice of the coordinates. The rays of an 
ordinary system being perpendicular to the surfaces which have for equation 


V = const., 

compose in general two series of rectangular developable pencils, and are tangents to two caustic 
surfaces. Let us therefore denote by %,, y,, z,, three rectangular coordinates so chosen that the 
axis of z, coincides with some given ray, and that the planes of a,z, and y,z, are the tangent 
planes of the two developable pencils to which that ray belongs; and let a, 8, y, denote, for any 
proposed ray of the system, the cosines of the angles which the ray makes with the axes of 
%,, Y 2. The equations (A) (B) (C) (D) (E) will apply to the coordinates thus chosen, by simply 
changing 2, y, 2, to ,, y,, 2,; and by changing y to its value 


EE 24 82 
y=V1—-@— ğ=1-Ž te — y” — y® — &c., 


in which 

1.3.5... (2i +1) (a?+ 8y 
TERE RETT Ce ee 

* [This is the general integral of (A).] 

t [If through the point x, y, z on the given ray we draw a plane perpendicular to the z-axis, cutting the latter 
at A,and through A draw a plane perpendicular to the ray, cutting it at B, then U= — p CB, where C is the point 
where the ray cuts V=0. The function S (p. 268) is an extension of this function U, including the initial as 
well as the final ray of an optical system. The function U was employed by Lord Rayleigh in discussing the five 
aberrations of von Seidel by the Hamiltonian method (Phil. Mag. 15 (1908), p. 681).] 
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the function W will in general admit of being thus developed, 
W = pW +5 (Ac? + BB?) + pW + pW + &e,, (F) 


W, A, B, being constants, and W®, W®, W, being rational homogeneous functions of the 
two small variables a, 8, of the dimensions 3, 4, ù, respectively. The constants A, B, are here the 
distances upon the ray, from the points in which it touches the two caustic surfaces, to the origin 
of the coordinates ,, y,, z,; and the terms proportional to a, 8, «8, disappear from the develop- 
ment of W, by the choice which we have made of these coordinates, and by the principles of the 
former Supplement.* In this manner the function U becomes 


U= pW — ps, +5 (2, +A) a? + (z, +B) B} +pU% + pU" + &e, (G) 


in which 
U i+ = Wears, Yit — Waits + gyei ; 


and the two first of the equations (D) become 
d d 
ama +(e Ay 9h; B=B, + (e+ B S (Œ) 
if we put for abridgment 


@ y 
at TEN B= 7B p=—(U% + UM + &.). (I) 
o $ do jet i 
n account of the smallness of — , 1p” the quantities a,, 8,, are approximate values of 


a, 8; and to develope a, B, HRT a or any function of them, F («, 8), in a series of ascending 
powers of these approximate values, we have, by the theorems of Laplace before referred to, 


F(t, B) = F, + ayo ig ia (i (z£) ie) asp (z (ie) ) (K) 


PG Fapa t [n+ TP Ce, + By 


maa) (32) +1 = (sey "E ($ E. Ta EIME- d (ay i 
n4 1p [n F pG, + Ay (2,4 Bye 


the functions F,, ¢,, being formed from F, $, by changing «, £, to a,, B,, and [n + 1]"*4, [n’4+1]}”’+1, - 
being known factorial symbols;+ we have therefore, 


drin 


, 


ona.+ gg dae (i) dard (daag te) (a8) ) 
ae Gray aa Cea cra ace aa 
aaz: (a8) ranr (deag (aa) (a8) ) 


E T .° 2,n da” dh,” 
B= B, (mo — T+ f += (n, n')0, 0 [n+] ntin z, + APH (z, + era” 


* (Cf. (D”), p. 127; the focal distances AERAN to ¢=0 and $=5- Hence, denoting partial derivatives 


by subscripts, W,,=0, R= — WaalVaa, Ro= — Weg/vgg. In these expressions W is a homogeneous function of 
a, B, y; but, since a=8=O for the given ray, these partial derivatives have the same values as when they are 
evaluated for W as a function of a, 8. The linear terms in the expansion vanish by (C).] 

t [See Appendix, Note 4, p. 468. ] 
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Now, if we differentiate V as a function of the three independent variables a,, 8,, z,, we have 
by (B) and (1), 


dV dV dV 
dq eT oP ap, ~ bP (+ 8), da ot Oe + BB, +7); (M) 
we have also V = wz, — pW, when œ, B,, vanish; and therefore, 


V= pz, — pW + p | ((¢, +A) ada, + (2, +B) Rdg), (N) 


z, being considered as constant in the integration,* and the integral being so determined as to 
vanish with œ, 8,. Substituting in this expression (N), the developments of «æ, 8, and perform- 


ing the integration, we find the following development for r ; 


r =z, — WO +4 {(z,+ A) a+ (e, + B) B) 


a ($ n+2 d” dd, n+2 
taS et ee a A „a a) 
a AAS RASTE EA OY ahd (E AEN 


driw ( do, (E Bi ($ yS 
43 00 da,” dB,” \da,dB, \da, dp, à (0) 


which is another form for the integral of the partial differential equation (A), obtained from the 
elimination (D). And if we wish to introduce any other rectangular coordinates æ, y, z, into the 
expression of this integral (O), instead of w,, y,, z,, we may do so by the known methods, by putting 
æ, = (æ — æ) cos. ææ, + (Y — Y,,) COS. yx, + (Z—2,,) COS. Z2, 
y, = (x —&,,) cos. xy, + (Y — Y) COS. YY, + (2 — Z,) COS. ZY, (P) 
z, = (x — x) cos. x2, + (y — Y,,) COS. yz, + (Z — 2,,) COS. 22, 
Ly, Yp» Zp, being the values of x, y, z, that belong to the point upon the ray which had been taken 
for origin. 


Verifications of the foregoing Developments. 


4. We may verify the form (O) which we have thus found for the integral of (A), by the 
following condition, resulting from (M), 


LIE AE. OE ON r = ee 

dep grAda p54 Bap, we e @ 
of which each member is an expression for the cosine y of the small angle which a near ray 
makes with the ray that we have taken for the axis of z,. The condition (Q) may be put under 
the form 


d V 


da, a Ce + BB) =VI— at — e, (R) 


* [The expression in the brackets {}is dV /p for a displacement in the plane z,=const.] 
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in which, by (O), 
d V a? +B? _ do, 
dọ, n+l do, d” dd, n+l dd, 
aun)" Or 
de i n+2 d” dọ, n+2 
nE AA aap (as) 
Zoo FOTODE Fay POTD) FBP 


qntn'+1 d?o, dg, n+l dọ, n’+1 
+ > w, 00 da,” dp,” ‘dz, Par dB, \da (Z ') (25) ) 
(mn) 0, OFT Ain +1)" G, + Ay (4, + By 


-sanit (Cot a sem ar (det (aa) (8) 
(n, wia ( n). [n] (2, + A)” (2, + By * 


d” dd, n+l d”? dd, n+l 
Be- aaa (38) 
w da; \da dp, 
arenas eat 
a dd, (d"d 
E doar te (2) GY”) 
+, Zen, n) 0, 0 njin + 1)" "3 (z, + A) (2, + By + 


, dawdae (itp (ae) (38) ) ; 


(n, n’) 0, omt iA mn e, FAs, F. aiin P e, FA)" G, + By 


(S) 


and, by (L), 


+B,> (T) 


while the development of 
VI-a- e 
may be deduced from the general formula (K) by changing 
F(a, 8) to y=VJ1 ate OP ome FFs 


F, to y=VJ1 —a,'— B,. 
To compare these several developments, and to examine whether they satisfy the condition 
(R), we are to observe, that from the nature of the function $, we have by the foregoing number, 


=- (yi ty ke) = eS em ter j (U) 
Pd, TE he Pd, =p 4%. C$, Oy, , 
da,dz, '"" dæ? dB,dz, “'' dB,’ da, dß,dz, = da, dB,’ 
and aan 


ina" @)-2. "a" 
n’+1 A a: “4 
+ (Bt) yA" + (00) CS GE 


and 
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by which means the difference of the developments (S) and (T) becomes 


a, n+1 Ki dp,\"** dy, 


da,” dB,” (dæ, dp, dp, dp,/__ 48," dB, \ da, dB 


, „o ta apr AED, dnu LARR aa ad a E px 
(n, n’) 0, 0 [n+ 1H [n T 1]” yi (z, + A)*t1 (z, + BY +1 ? 


(V) 
and the series in this second member being exactly that which would result as the development of 
y=VJ1-@— f, 


from the formula (K), we see that the condition (Q) or (R) is satisfied, and the sought verification 
is obtained. 


Another verification of the foregoing developments may be obtained by applying the general 
expression in series (K), for any function F of the cosines a, 8, to the case where this function 


18 = i . We find, first 
da 


1 db ys nf (ay ae) isp (Gey set) | 


43 „o, ada hap ldaždß,\da,) _\JB, daa \d8) dB, \da,) t dadp, da,’ \d8 
(n,n) 0, 0 [n+ 1]? [n? HLH (e, + A) (z, + By 
(W) 
which may be put under the form 
aS jana (tarda, ae) (as) ) 
dp y 2 dam \da) y „<o dandB” \da dB, \da,/ \dB, X) 
da ° [nt Ae FAS OA % nn ETP Cz, + Ay (2, By 
that is, by (L), 
$ = (¢,4.A) (@-a), Y) 


which agrees with the conditions (H). A similar verification may be obtained from the same 
conditions (H), by considering the development of n : 


Finally, we may observe that the condition 


p= + By tya- =a, + By ~ 2 (Z) 

becomes, by (G) and (I), 
ys a B: z r 
we - Wo +(e, + A) (aa, 5) +(+ B) (88, -Ẹ) +; (A’) 
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drin a ad, ($ pe (58 y B ad, (se y +1 d (ey d*¢, ($ i d ($ a 


> 
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in which, by (K) and (L), 


n+1 d” (E ka driw Paa dd, (= p. ($ y H) 
a ee = ig ah. n +2 da” \da S », œ da,"dp” da,dß,\da dp, 
a (n)0 n+] + I era (n, n’) 0, 0 [nin + 1)" AG, +A (z, + By , 
n+l da" ae dew 7 dẹ, (dpm ig) 
BB ae i gr _y 2 n+2de," (73 _y 3 dæ” dB,” us na (2 ia) (a3, 


a (2) ($ jv 
$ Pa $ + = da,” da, dp,” dp, 
/ (n) 0 [n+ PG, +A)’ 1} (z, +A (n+ IPF ip (2, + By ti 
drn do, dọ, n+l do, n’+l 
Aw +3 ganda (dada, (dar) (as) ) sen 
[n + 1p" [nV +1)" +3 (2, +A)" (z, + By R 

so that we are conducted by this other method to the same expression (O) for the characteristic 
function of an ordinary optical system, as that which we before obtained by performing the 
integrations (N). In all these expressions the sign =,,, ,,.%y ¢ denotes a summation with reference 
to the variable integers n, n’, from zero to infinity. 


+ Benny 0,0 


Case of a Plane System. 
5. A similar analysis may be applied to integrate the partial differential equation 


dV\?  /(dV\* , 
(ae) + (Ge) = (©) 
to which the equation (A) of this Supplement reduces itself, when we consider a system of rays 
of ordinary light, contained in the plane of xz. In this case, if we put 
L, = (& — @,) COS, wu, + (Z —Z,,) COS. 2%, 5 
(D’) 
2, =(@ — 8) COS. %2, + (2 — Z) COS. 22,, 


we may suppose «,, 2, to be new rectangular coordinates, in the same plane as g, z, and such that 
the axis of z, coincides with the direction of some given ray of the system: and we may denote 
by @, y, the cosines of the angles which any near ray makes with these new axes, so that 

y= y 1- a. 
We shall then have for one form of the integral of the partial differential equation (C’), the 
following combination of equations : 


av _ av 
pa, = Ts V+U=a—, (B’) 
between which g is conceived to be eliminated, and in which 
2 
U= W —'pyz, =p W0 = ps, + PEDE u; 
-=E go (WO + 2,949) ; (F’) 


1.3.5... (2+1) a*t. 
2.4.6... (+2) +A’ 


HMP 20 


WHI = gtt3 wis; yat = 
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W, wis, being constant coefficients in the development of the function W, according to the 
powers of a, and A being another constant in that development, namely, the distance upon the 
given ray, from the point where it touches the caustic curve of the plane system, to the origin 
of æ, and z,, The first of the two equations (E’) becomes 


when we put 


a= Z, + z.+A 3 
fen gives therefcre, by the well-known theorem of Lagrange, for functions of a single variable, 
d” dF, (dd, n+l 
DEAA 

mo mI PAE tI 
F(a) denoting any function of a, which admits of being developed according to positive integer 
powers of a,, and F,, ġ,, being formed from F, ¢, by changing a pa, a,. The cosines @, y, may 
therefore be thus developed, 


F(a)=F,+> (G’) 


d” (3 a 
» da” \da, 
«= 4 +200 FIG rA ge 
a dy, (y 
L 43 2 da da, dey 
TSU thw Ce TPs, + APH 
if we put 
y=V1- a’. 
And since V may be thus expressed, 
V=pe,—pW + p(z, +4) |" ada, q’) 
because 
dY = p (ada, + ydz,) = pa (z, + A) da, + u («, +y) dz, 
and because V becomes uz, — pW% when g, =0, we find, finally, 
d” (2 fa 
aie (2,+ =le a, da,” da, š 
ea o e Rs Ta? tae 
H = Pony +9, + Zin) 0 [n+ IP tay Z+A yur (K’) 


This form (K’) for the integral of the partial differential equation (C’), may be verified by 
observing that it satisfies the condition 


t ay E) 
V being differentiated for z, and a, as two th ee variables ; because 
a is +D) as = (tees. 

gana (Fe) agea (E ted gee (GB, cal’) 
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and therefore, differentiating (K’) as if a, were constant, 
d” ($ hy d” dy, (= n+1 
Pave as du \da o ap ae ae) 
p dg, ute, +E 26 mre rA t wo n IP G F Ay 
that is, by (H’), 


N’) 


Case of a System of Revolution. 


6. Another particular case of the partial differential equation (A) deserves to be considered 
specially; namely the case of systems of revolution, symmetric about some single ray. In this 
case, if we take for the axis of z, the ray which is the axis of the system, V will be a function of 
z and of a + y?; and if we put 


a +y’ =n, (0°) 
we may in general suppose V developed according to positive integer powers of n, in a series 
satisfying the condition, 

dV dV s , 
(a) te (a) iia ®) 


To integrate this partial differential equation (P’), which is a particular case of (A), we 
may employ the corresponding case of the general system of equations (D), (E), putting for 


abridgment 
. @ +A =e, (Q’) 


and considering the quantities W, U, as functions of e, which we shall suppose capable of being 
or according to positive integer powers of that variable. In this manner we shall obtain 


dU dU dU 


da de? dp = 28, | nr 
and therefore by (D), an 
2 
pin = te (T) V+U= 27. (S’) 
We have also by (E), 
U=W-ypzv1-—e, oe Get (T’) 
in which ee 
; y=V1—e; 
and we may put 
Ww wo pAe +2 ° 
=p +> tHo“ Wei+a 5 
dW pA 
ee + yo (G+ 2) eft Waite; 
(U’) 
dU zt+A o i 
PA a )— pġ; -p= go Uime; 


1.3.5... (2i+1) 2. 
Tnm 64+ 2) ware 3 a6... (0422 
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in which W, wzi+4, A, are constants of the same kind as before, A denoting the distance of the 
origin from the focus of central rays. Hence, if we put for abridgment, 


“(+A (@+A/y’ id 
e, is an approximate value of e, and we have* the following relation between e and ¢,, 


dey herp” , 
z+A GFAP’ w 


which gives, by the theorems before referred to, 


nti df, ear? \ "4 
P TE Sess milihna 


e=e€,+ 


(X’) 


J(e) being a function of e, and f,, ,, being formed from f, y, by changing eto e. We have also, 


by (S’) (T”) (0’), 
Fac Woi f(e e Ta) de 
m H Jo \de de? 
=e- WO EEA) dens f yde, S 
0 
and therefore, by (X’), 
V zt+A ' €, 
p =Z -WO + Cnet. 2er, + | yde, 


o 4n+1 d” Z+ A dy, ( er? say" 7 , 
ii Zeno [n + 1p (z Fi Ayt de,” (a a 2e, de, ri Y) e — z+ A (Z ) 
in which, 
=—2 5 tty. jet SPs i. 
Vi =— 20 Unie; de, “0 (i+ 1) Uime; 
€, an 7] 
f e T 553 “= e+ AP 


The development (Z’) is one form of the integral of the partial differential equation (P’); 
another form of the same integral may be obtained fron the expression (K’) for the characteristic 


function of a plane system, by changing a, to Ag A? and z, to z, and supposing w;43 = 0, and is, 
V a EFASE a 1+5 do \”+2 
-=z Wo 43° o ( (=) ), ye 
pO Wa ath as eee Gaal’ Gay el 
in which 
o E ý 1 Beo ME A OA E i D UET B” 
es -Zos gy ayen (Vn 2.4.6.. I) vo 


Each of these forms gives, when we neglect 74, the following approximate expression for the 
characteristic function V of a system of ordinary rays, symmetric about the axis of z, 


=z- wor 7 —_ (2 + Bua) 9? (2+ 16w)n? | (2 + Bw)? rè, 3 
aT + 2+) S(z+A)* 16(2+A) 8(z+ Ay’ (C 


* [By (S).] 
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in which, y= a+ y*, and W®, A, w4, we, are constants in the development of the connected 
function W, such that when we neglect the eighth power of the sine of the angle contained 
between a near ray, and the axis of revolution of the system, we have 


A = WO + eer’ + Ws (a? + B$ + We (a? ihe B$, (D”) 


a, 8, being, as before, the cosines of the angles that the near ray makes with the axes of æ and y, 
to which it is nearly perpendicular. 


Verification of the Approximate Integral for Systems of Revolution. 


7. The approximate expression (C’’) for the characteristic function of an optical system of 
revolution, admits of extensive applications: it is therefore useful to consider other methods, by 
which it may be obtained or verified. An immediate verification may be derived from the partial 
differential equation (P’) of which (C’”) ought to be an approximate integral; namely, by com- 
puting from (C’’) the approximate expressions of 


1 /dV\? 4m (dV 

p’ (a) * omens it (an)? 
and trying whether their sum is unity, when 7‘ is neglected. Putting for this purpose the 
expression (C’’) under the form 


ae ai cg a Ta 

pat ~ WO +oG ea) Sa+APt ios ap 
(A — 8w) n? ; (32wa — 16we — 3A) n? , (A — 8wa)? n? Œ") 
“8+ Aye 16 (2+ 4) SkA ' 


we find, by differentiation, 
EEF ei 019 a ANORA.. h 


pd” XCEW. SA 8(z+ AÝ 16(2+A)® 
(A — 8w) n? 3 (32w, — 16w — 3A) n? 7 (A — 8u4)? 9°. 


~ 2+ AP 8 (2+ A) 8+ AP? E) 
2dV _ 1 Pade, dubs! 3n? (A — 8w) n 
ud z+A 2@+Al*8G+Al* 204A 
3 (82w4 — 16w — 3A) 7? | 3(A — 8x4)? 9? , 
8 (z+ A)? 4(z+Ayr 
and therefore, neglecting 7‘, 
4(%)- SOL, aN ys. ANS IE Frat ha 
p? A (2+ A)? t+ GHA + A¥ 
ee (11A — 112w, + 48w) n? 7 (A - 8w,’ n? 3 
i 4(z+ A)? ~ 4&@+Ap ’ (a) 
4n (7) - mm, (A-8u) 9? 
p? \dn GFAP “+AT (z+A) (2+ AP 


ge (11A — 112w,+ 48wW¢) n? 7 (A - 8w) n? S 
4 (2+ A) + 4G@+ApPp ’ 
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expressions of which the sum is unity, as it ought to be. We may remark that the former of 
these two expressions represents the square of the cosine, and the latter the square of the sine, 
of the angle which a near ray makes with the axis of revolution of the system. 


Other Method of obtaining the Approximate Integral. 


8. Again, the approximate integral (C’’) of the partial differential equation (P’), may be ob- 
tained in the following manner. Since V is supposed capable of being developed according to 
positive integer powers of n, let us assume 


anes Vo +V n + ya n? + yea 1, (H”) 
H 


neglecting 74, and considering V ®, V®, V®, V ®, as functions of z, of which the forms are to be 
determined. To determine these forms, we have, when » = 0, 


dV os a 
dV ava @vy avo rT dV ei ETA 
dz dz’ dnde H de > dpde V de > dade Pe ’ 


(r) 


The equation (P’) shews that a + u, when y =0; and =A is positive, if we suppose the 


motion of the light directed from the negative towards the positive part of the axis of z; we have 
therefore, by (I’’), 
dV T 
qpa i. (K’’) 


The equation (P’) gives also, by differentiating it with respect to », 
_dV BV dV d?V (G): 


T o hra r E omens 
caret yey AVV (BV dV av. i 
-T artt (a) +a at tl aes) +8 ay ah? od 


=W OV o OV OV oy AVY gt ee ht +12 (Z5); 
= dz dde |” ddz dnde ©" dn dnt" da? dy” dy dnp dn? 


and, making y = 0, we find by (I”) the following equations in ordinary differentials, from which 
V%, V®, V®, are to be deduced : 


o=% Lavan, 
de 
a) 
o=% +syo ya +4 (Are, M”) 
ave dvuave 
SS eee a 3) 4 (2)2 
o=% +1270 yoo 2" saya, 
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These three differential equations, when divided respectively by V %2, V%4, V8, can easily 
be integrated, and give, when combined with the integral of (K’”’), 


Vo = VO—z; 
1 
ETERNA N”) 
ya 
a= aa t 22; 
Aue 4 (v + V2)", 
yo Sem 


Vo, V1, Va, Us, being the four arbitrary constants introduced by the four integrations. The functions 
Vo, V®, V®, V®, are therefore of the form 


1 Ug — 2z 
 — . 2 EVE o CAR po R a 
V PF, V ztu’ V (zF 4)" 
4z +u | 4(v, +02). 
(22+%)® © (2¢+%)? ’ 
and these forms for the coefficients of the development (H’’), agree perfectly with the development 
(C’”) or (E”), when we establish the following relations between the constants : 
w=- W0; wn =24; v= — 164; vs = 8 (16w, — 8w — A): (81) 
we see, therefore, that the present method of integration confirms the former results. 


(0”) 
V3) = 


Connexion of the Longitudinal Aberration, in a System of Revolution, 
with the Development of the Characteristic Function V. 

9. To give now an example of the optical use of the development which has been thus 
obtained, let us consider its connexion with the aberrations of the near rays, from the principal 
or central focus. We have already remarked that the constant A denotes the distance of the 
origin of coordinates, upon the central ray, beyond this principal focus, in such a manner that 
the focal ordinate is =— A. For the ordinate Z, of intersection of any near ray with the central 
ray, we have by the fourth of the equations (C), of the present Supplement, 


1 èw ” 
7= (+B) by’ me? 


if we form the coefficient a by putting W under the form of a homogeneous function of a, £, y, 
of the dimension zero, with the help of the relation a? + 8? + °= 1, and by then differentiating 


this function, as if a, 8, were constant, and y the only variable. Employing therefore for x 
the development 


w A (+B) y e Wara (a + Bt £ 
MOF See BR oft OA (R) 


which is of the homogeneous form required, and, after differentiating for y, making a + 8? + y’ =1, 
we find for the ordinate Z, 


Z=-A+A (Ly) -Ypo (2i + 4) ws (@? + BM, (S”) 
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a series of which the term — A being the ordinate of the central focus, the remainder is the 
longitudinal aberration: y is the cosine of the angle which the near ray makes with the central 
ray, and a? + £? is the square of the sine of that angle. If therefore we denote the aberration 
Z + A by A, we may develope A in a series of the form 


A = L (a + £?) +L, (a? + 6)? + &e., (E) 
L=}4— 4u, [,=44 +2w,— 6w. (U”) 


And if by these relations (U’’), we eliminate w4, we from the approximate expression (E”), 
we find the following formula : 


0 y L 
V=n@-W) +57 ay SEN cee tg 3) 


me. aabb. Oe h 
CoH z+A + erupt it ip 


which shews the connexion in a system of revolution between the development of the longitudinal 
aberration A, and that of the characteristic function V. 


in which 


+ 


Changes of a System of Revolution, produced by Ordinary Refraction. 


10. Suppose now that the rays of this system of revolution fall upon a refracting surface of 

revolution, having for axis the axis of the system, and having for equation 
Z = Zo + 21N + 220? + 257° + &e,, (W”) 

in which 7 is still = æ? + y? = the square of the perpendicular distance of a point æ, y, z, from the 
axis; and let yw’ be the refracting index of the new medium into which the rays pass after 
refraction. It is evident that in this new medium, the rays will compose a new system of revolution, 
symmetric about the same axis as before; and we may in general suppose the characteristic 
function V” of this new system, which is analogous to V of the old, developed in a series similar 
to (V”), 


9 L’ 
_wW'o AEE ss, tend ec ~ ) 
AP RA Oto Hint 4 (2+ A’ 4+’ a ie 
pe’? (20, —5L’ 6L” 4 
tiera zea Grp I ©) 
the constants A’, L’, Li being similar to A, L, Z,, in such a manner that the ordinate Z’ of 
intersection of the axis with a near ray, is 
Z'=— A'4+L' (a? + B’) + Li (a? +B), (x) 
if a’ + 8” denote the square of the sine of the angle which the near ray makes with the axis, 
and if we neglect the sixth power of this sine. To connect the new and old constants in the 
development of the characteristic function, we have, by the nature of this function, and by the 
principles of my former memoirs,* the condition 


0=AV=V'_-YD; (Z’’) 


which is to be satisfied for all the points of the refracting surface, and which may therefore be 
differentiated, considering AV as a function of z, n, namely the difference of the developments 


* [See p. 117.] 
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(V’’) (X”), and considering z as itself a function of 7, assigned by the equation of the refracting 
surface (W’’). In this manner we find, transposing the symbols A, d, 
dV dz,dV 
0= "AT +7 di A Tz’ 
M ds, @V aNs Ee ds aY. 
( o 
noy ds , ®V ds\! @V . (ds\?, &V 
ae al “Adlets dedm +3(F =) ata; + (ae) aap 
ds, @V dedy Pri dz =} 
and, making after the differentiations » =0, we have 


Z = 20; ira, 9. TAi on; os = 62s; 
dn dn? dn? 


dV , CV: BV i 
ail Ao5=0; AT5=0; 
dV ae ae PEND jaan m 
-$A sn =}A GFA Adam’ GA (RD 
arali ESEE f A av a {= ape Se 3u 
IGA FGE+AP? Oded O GFAP 44 Ay? 
Atv m 2L, — 5L 3I 
aap TA araplacea* Gray a): 


We have therefore, Ist, for the change of — A, the ordinate of the central focus, 


Om Aa + 2a du: i (0) 
IInd, for the change of L, the first or principal coefficient of aberration, 
on 2uL land Ép t A SLP > r 
OAS Ay GAP A er ap t An: i, 


IIIrd, for the change of Lı, the coefficient of the fourth power of the sine of the angular aberration, 
in the expression of the longitudinal, 


a fetter. “am 2L 
ny a ofa (stea*aray*!)* aa aati) 


+322 A —32,A. +62A4p;  (E”) 


‘Gray Gray 
z being the ordinate of the point of central incidence. With respect to the present meaning of 
the sign A, we may remark, that the first of the three equations (C’”’) (D) (E’”) is equivalent 
to the following: 


FEI +2 = ypt Qaim’; (F’”) 
and the two others are to be similarly interpreted. 


HMP : 2I 
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Example; Spheric Refraction; Mr. Herschel’s Formula 
for the Aberration of a thin Lens. 


11. These general equations for refracting surfaces of revolution may be adapted to the case 
of a refracting spheric surface, by making 


1 u 1 ttt 
43> or? a= 33 23 = 7675? (G ) 

the two first, for example, becoming * 
‘ued p H - mt 
o=a(-* +4); (B) 


O=A (4 H 2u 5) : 

(¢+A)* (¢+AP® r(z+AP r’ 

which contain under a convenient form, the known theorems for the change of the central focus, 

and of the principal coefficient of aberration, by refraction of a spheric surface; r being the radius 

of this surface, and being considered as positive or negative, according as the convexity or con- 
cavity is turned towards the incident rays. 

If, for instance, we consider an infinitely thin lens in vacuo, having yp for its refractive index, 
and having r,r’, for the radii of its two spheric surfaces, (positive when those surfaces are convex 
towards the incident rays,) we may take the point of central incidence for origin, and the equation 
(H) will become, 


a”) 


Ae wine h MEL bs 1 SE 2.1 Tp tr 
Ora Sat ae perce coe 
— Á, — A’, — A”, being the ordinates of the central focus in the three successive states of the 
system; and similarly, (I’’’) will give 


o= HL! a -(4 3) -( 2u a) tS 


Bie tate. E ee rā rA? iai 


o= 22 _ 2E ae ] H 2 2u l=p. 
= gi — “ae ~ (gra fa) T O + Ga 
L, L', L”, being the three successive values of the principal coefficient of aberration. Adding the 


two equations (L’’’), the intermediate coefficient Z’ disappears, and we find, 


1) TE VE Meh I qs wae ie | i 
0= ga-nt) eaaa) t e-h a-a) OL") 
in which, by (K’”), 


(Eid) 


1 fe 1 p—l. 1 AN 1 1 > rrr 

aa es m a e- Dla) yar, 
and therefore, 

Pn Mia daak, aT M® ye j 

ga E py) re e eet 


* [The quantity in the brackets in (H”) is the so-called “Zero-Invariant” of Abbe; cf. J. P. C. Southall, 
Geometrical Optics (1913), p. 159, where the introduction of the A notation in the theory of aberrations is 
incorrectly attributed to Abbe.] 
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if we put for abridgment 


- 2p + u? 
14+2u—-2y2 p 
a . 
PERA i T ANTATT] 

443- ene (P) 
Moe tse - 


M® =24+3. 
p 


It is easy to see that the formula (O’”) agrees with the expression for the spherical aberration 
of an infinitely thin lens, which Mr. Herschel has deduced by reasonings of a different kind, in his 
Memoir “On Aberrations of Compound Lenses and Object-glasses,” published in the second part 
of the Philosophical Transactions for the year 1821; and in his excellent “Treatise on Light,” 
published in the Encyclopedia Metropolitana.* 

The elegance of this formula of Mr. Herschel, and the important consequences which he has 
obtained from it, have induced us to shew how the same expression may be derived from the 
development of the characteristic function of an ordinary system of revolution, assigned in the 
present Supplement. The same form of development, and those other forms which we have assigned 
in the same Supplement, for systems not of revolution, contain the solution of other optical 
problems, of which we hope to treat hereafter.+ 


* [Vol. tv, pp. 341-586. Although the title page bears the date 1845, the Encyclopædia appeared in parts 
during the years 1829-1836. Herschel’s article is dated Dec. 12th, 1827, and mentions Hamilton’s “Theory of 
Systems of Rays,” which was being printed at that time.] 

+ [See footnote to p. 165.] 
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